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ABSTRACT 

We demonstrate that a natural consequence of an asymmetric kick imparted to neutron stars at birth 
is that the majority of double neutron star binaries should possess highly eccentric orbits. This leads 
to greatly accelerated orbital decay, due to the enormous increase in the emission of gravitational 
radiation at periastron as originally demonstrated by Peters (1964). A uniform distribution of kick 
velocities constrained to the orbital plane would result in ^24% of surviving binaries coalescing at least 
10,000 times faster than an unperturbed circular system. Even if the planar kick constraint is lifted, 
^6% of bound systems still coalesce this rapidly. In a non-negligible fraction of cases it may even be 
possible that the system could coalesce within 10 years of the final supernova, resulting in what we 
might term a "double supernova" . For systems resembling the progenitor of PSR J0737-3039A, this 
number is as high as ~9% (in the planar kick model). Whether the kick velocity distribution extends 
to the range required to achieve this is still unclear. We do know that the observed population of 
binary pulsars has a deficit of highly eccentric systems at small orbital periods. In contrast, the long- 
period systems favour large eccentricities, as expected. We argue that this is because the short-period 
highly eccentric systems have already coalesced and are thus selected against by pulsar surveys. This 
effect needs to be taken into account when using the scale-factor method to estimate the coalescence 
rate of double neutron star binaries. We therefore assert that the coalesence rate of such binaries is 
underestimated by a factor of several. 

Subject headings: pulsars :binary, gravitational radiation, merger rates 



1. INTRODUCTION 

Beginning with the detection of the Hul se & Taylor 
double neutron star l|Hulse fc Tavlorl 119751) . pulsar as- 
tronomers have unveiled a small but growing sample of 
double neutron star binaries over the last thirty years. 
Owing to the extreme regularity of their radio pulse 
emissions, these exotic objects have proven to be in- 
valuable scientific laboratories, particularly for tests of 
general relativity, as they are some of the few known ob- 
jects for which the effects of gravitational radiation are 
measurable. In fact, the first confirmation that gravi- 
tational radiation existed was through measurements of 
the orbital decay of the fi rst known double neutron star 
ijTavlor fc WeisherdlT989T) . 

Recent surveys have started to unveil a better picture 
of the double neutron star population. The Parkes multi- 
beam surveys have discovered three probable double neu- 
tron star systems llLvne et al.ll2000t iBurgav et al.ll2003l 
iFaulkner et al.l2005fl. and Arecibo surveys have det ected 
two more (Wolszczan 1991; Champi on et al.ll2004|) . Al- 
though we are still dealing with small number statistics, 
the population is now large enough to see trends emerg- 
ing. 

A phenomenon of crucial importance to the binary co- 
alescence rate is the kick velocity imparted to a neutron 
star during its supernova, a dominant factor in deter- 
mining the eccentricity and coalescence time of the sub- 
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seque nt orbit l|Hillslll983t iDewev fc CordesHl987t iBailesl 
1989). Binary neutron star formation and evolution also 
depends greatly upon this little-understood supernova 
kick, since a double neutron star system must survive two 
such kicks in order to remain bound. To some extent, 
knowledge of the kick velocity distribution can be ex- 
tracted statistically from measure ments of the observed 
samp le of double neutron stars ( Wille ms fc Kalogeral 
2004). The magnitude of the kicks can be inferred by 
the space velocities of pulsars, which are an order of 
magnitude higher than that of their progenitors, the OB 
stars. Investigations of the pulsar population demon- 
strat e that the averag e kick must be of order 400-500 km 
s" 1 llLvne fc Lorimerill99l: ICordes fc Chernoflll998|) . 

Further evidence for pulsar kicks comes from preces- 
sion e ffects observed in the binary pulsar PSR J0045- 
7319 (Kasoi et al. 1996) and spin-orbit coupling in the 
relativistic binaries PSR B1913+16 PSR B1534+12 
and PSR ,11141-6545 llWeisberg et aTM989l iStairs et alJ 
l2004tlHotan et alJ2004|) . Assuming tidal forces and mass 
transfer align the spin and angular momentum axes prior 
to the ultimate supernova explosion, these systems sug- 
gest a kick out of the orbital plane has occurred. 

Binary pulsars are in some sense fossil records of the 
progenitor system from which they evolved. We date 
these systems by the characteristic age of the pulsar, and 
can therefore extrapolate back to the orbital conditions 
immediately after the final supernova explosion. What 
prohibits us from going further back in time to derive a 
unique progenitor is the possibility that the newly-born 
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pulsar received an impulsive kick at birth. Since the kick 
velocity distribution is largely unknown, including its ori- 
entation with respect to any reference point, such as the 
pulsar spin axis, we are left with a rather unsatisfactory 
family of potential progenitors for any observed neutron 
star binary. 

The relativistic nature of these objects has brought 
attention to the galactic rate of double neutron star co- 
alescence - the merger event at the climax of a long or- 
bital decay and inspiral due to gravitational radiation. 
Such events are some of the most significant anticipated 
sources of gravitational waves, detectable by interferom- 
eters such as the L aser Interferometer Grav ity- Wave Ob- 
servatory (LIGO) ijAbramovici et all :19921 . Double neu- 
tron star coalescence is also considered a major candidate 
for explaining the origin of a subset of the mysterious cos- 
mic gamma-ray bursts observed by gamma-ray observa- 
tories such as the BATSE instrument abo ard the Comp- 
ton G amma-Ray Observatory (CGRO) l|Meeean et al.1 
1992). 

Much research has been conducted into pre- 
dicting the double neutron star coalescence rate, 
using Monte C arlo p opulation synthesis models 
(iBelczvnski et all l2002t iTutukov fe Yungelso"nl |1993t 
iLipunov et alJ Il997t iPortegies Zwart fe Spreeuwl Il996fl 
or by statistical extrapo l ation fr om the observe d sam- 
ple llCurran fc Lorime^l 119951 iKalogera et"aT] 120011: 
iKalogera et alJ l2004|) . Monte Carlo methods are cur- 
rently hampered by the vastness of the population 
parameter space and compounded uncertainties, and 
statistical methods suffer from the small sample size of 
known double neutron stars. 

It is customary to estimate the event rate that ad- 
vanced LIGO will see by estimating what fraction of 
neutron star binaries we could have detected in pulsar 
surveys if they were of similar period and luminosity to 
those known, and assigning a "scale factor" to each. We 
then multiply each member of the population by their 
scale factors and obtain an estimate of the number of 
similar systems in the Galaxy, thus obtaining a coales- 
cence rate for the Milky Way that can be used to infer 
an event rate for the local Universe by estimating the 
number of Galaxies out to the LIGO detection limit. We 
will refer to this method as the scale factor approach. 
Whilst limited by small number statistics, the growing 
population of observed neutron star binaries makes us 
increasingly confident that coalescence is relatively com- 
mon in the Universe. 

The only alternative approach is to attempt to simu- 
late the entire binary pulsar population. Uncertainties 
in initial conditions and the physics of binary evolution 
and kicks mean that such estimates vary by mor e than 
two orders of magnitude l)Belczvnski et al] l2002). The 
population synthesis method is therefore of limited use 
in estimating the merger rate. 

For now, the coalescence rate is better estimated by 
considering the scale factor method, however when this 
is done we find the coalescence rate is dominated by th e 
double pulsar PSR J0737-3039A ijBurgav et al.l f2003). 
This serves to demonstrate the heavy influence of small 
number statistics on the uncertainties in the scale factor 
method. 

The scale factor method also suffers from the extent to 
which the relativistic binary pulsar population remains 



invisible to pulsar surveys. For instance, if there existed 
a population of neutron stars that were completely invis- 
ible in the radio we would have no way to estimate their 
population from radio observations. Similarly, if only 1% 
of recycled neutron stars emitted radio waves, we should 
multiply the scale factors currently in use by 100. 

Our best estimates of the inspiral rate require us to 
think of all the reasons why double neutron stars might 
be invisible to our surveys, and take them into account. 
What we explore in this paper is a selection effect that 
removes binary pulsars from our surveys and has received 
little attention to date. This is the production of highly- 
eccentric short-period binaries with very short lifetimes, 
produced via favourably-oriented kicks in the binary pro- 
genitor. These systems are not long-lived, and the coa- 
lescence rate of the scale-factor approach needs to be cor- 
rected for it. Whilst this effect has been a lluded to from 
time to time (jTutukov fc Yungelsorfl 993) , we determine 
its importance for an arbitrary binary in the presence of a 
kick velocity that can occupy any of the phase space that 
will produce a bound orbit, as in this case the problem 
is scale invariant. 

In this paper, we derive the expected distributions of 
orbital eccentricity for any surviving binary pulsar, if the 
kick velocity is capable of uniformly populating the avail- 
able phase space after the explosion (Section Q). In Sec- 
tions and 01 we look at the lifetimes of the surviving 
systems due to gravitational radiation losses and iden- 
tify a population of binaries that coalesce very rapidly. 
These results are used to explain the observed distribu- 
tion of binary pulsar eccentricities (Section which at 
long orbital periods favour high eccentricities, and at low 
periods, small eccentricities. In Sectional we argue that 
the coalescence rate should be revised to take this into 
consideration, particularly for the double pulsar. 

2. ECCENTRICITIES 

In this section we examine the expected orbital eccen- 
tricities of relativistic double neutron stars from a purely 
analytical point of view, avoiding any assumptions about 
specifics such as the distributions of stellar masses, kick 
velocities and orbital separations. Instead, we present a 
study of the potential eccentricities of an arbitrary dou- 
ble neutron star system in units scaled to the circular 
orbital speed, such that the conclusions drawn here can 
be considered common to the evolution of all double neu- 
tron star systems. Indeed, all of our results also hold for 
any two compact objects (black holes, white dwarfs or 
neutron stars) where the final compact object (a black 
hole or a neutron star) is created in a supernova and 
receives a kick. 

We view the effect of kick velocity through a map 
of the second neutron star's post-supernova "velocity- 
space" - our representation of all possible relative veloc- 
ities which correspond to bound orbits. The boundary of 
this velocity-space region is circular, with radius equal to 
the escape speed of the system (y/2 in units of the circu- 
lar orbital speed) . Each point in the map corresponds to 
one possible orbital velocity vector, where the position 
(0,1) would correspond to the circular orbital velocity. 

To calculate the semi-major axis and eccentricity of an 
arbitrary binary, we derive formulae fro m two expre ssions 
for the total energy E of the system l|Hills! fl983'). and 
two expressions for the total angular momentum L of 
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Here /i = (Mm) / (M + m) is the binary system's re- 
duced mass, Mt its total mass, a its semi-major axis, 
e its eccentricity, v the magnitude of its relative orbital 
velocity oriented at an angle 9 to the plane defined by 
the vector that connects the two stars, and R the initial 
orbital separation. From these four equations we derive 
expressions for the semi-major axis and the eccentricity 
of an arbitrary binary, given the initial relative orbital 
velocity v: 



R 



(5) 



= [1 - v 2 cos 2 9(2 - w 2 )]^ 



(6) 

These formulae are applied here to generate a contour 
map of eccentricities in velocity-space. Note that the re- 
sults are independent of the choice of R. A code was 
written to uniformly populate the velocity-space region 
with grid points, then apply Eq. EJto calculate the eccen- 
tricity corresponding to each grid point (where each grid 
point represents one possible velocity) . The 2D data was 
finally plotted as a contour map, and is presented here 
in Fig. □ 

In interpreting this map, note that the relative veloc- 
ity vector just prior to supernova would be a little larger 
than 1.0 in the scaled units used here, because of the 
larger stellar mass prior to supernova. (For example, a 
pre-supernova mass of 2.0 Mq corresponds to an initial 
orbital speed of 1.102 units.) Also, for a sense of scale, 
note that the relative circular speed of PSR J0737-3039A 
is of order 600 km s -1 . We operate on the assumption 
that each system's final supernova is accompanied by an 
asymmetric supernova kick, which has the effect of dis- 
placing the initial relative velocity to a new position in 
velocity-space. This results in a new eccentricity and 
semi-major axis, which can be calculated using the new 
velocity and Eq. [S]and||Jl We further assume that the 
outgoing supernova shell's impact with the companion 
has negligible effect. 

Now, the most striking feature of the map in Fig. ^ is 
the very small size of the low eccentricity regions. This is 
surprising, as the currently observed relativistic binaries 
tend to have eccentricities below 0.3 and therefore must 
have been kicked into a very specific region of their post- 
supernova velocity-space (the third inner-most circle in 
Fig. This seems an unusually fortuitous coincidence, 
requiring a very consistently precise kick. Furthermore, 
with mean kick magnitudes of order 400-500 km s" 1 , we 
expect systems to be highly prone to being kicked out of 
these tiny low-eccentricity regions in velocity-space. 

The 2D eccentricity data generated for Fig. ^is pre- 
sented as a histogram and cumulative distribution in Fig. 




Radial velocity component 

Fig. 1. — Contour map of orbital eccentricity in the post- 
supernova velocity-space of an arbitrary binary star system, in 
which one star receives a randomly oriented kick. The velocities 
are measured in units of that required to place the post-supernova 
system in a circular orbit - for example, the relative circular speed 
of PSR J0737-3039A is of order 600 km s _1 . Prograde orbits have 
positive tangential velocities. The circular boundary of the region 
represents the escape velocity of the system, which is a/2 times the 
circular speed. 
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Isotropic kicks 




Orbital eccentricity, e 

Fig. 2. — Eccentricity distribution histograms are shown for bina- 
ries that are produced from randomly-oriented kicks which fill the 
available phase space for the planar (solid) and general (dashed) 
cases. In the general case the kicks can be out of the plane. Also 
shown is the probability that the final eccentricity will be less than 
a given value, which is simply the integral of the eccentricity dis- 
tributions. Note the high probability of eccentric systems resulting 
from a supernova kick in both cases. 



13 Here we have also provided eccentricity data for a 3D 
version of our planar velocity-space, uniformly populat- 
ing a "sphere" of radius \pl. This allows us to investigate 
the effect of supernova kicks which are not restricted to 
the orbital plane. 

The eccentricity histogram clearly illustrates a heavy 
bias towards very high eccentricities in the planar 
velocity-space, with over 50% of all eccentricities greater 
than 0.8. The distribution is less skewed when we allow 
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for kick velocities out of the orbital plane, but the bias 
is still towards very high eccentricities. We assert that 
this illustrates a strong inherent tendency towards highly 
eccentric systems arising as a consequence of supernova 
kicks. 

This finding is not reflected in the number of highly ec- 
centric relativistic binaries in the observed sample of dou- 
ble neutron stars, which tend to have low eccentricities. 
Analysing this from the standpoint of kick velocities, 
the observed systems seem to have been kicked consis- 
tently into a very specific region of their post-supernova 
velocity-space, the "bullseye" seen in Fig. ^ Our expla- 
nation of this tendency is discussed further in the follow- 
ing sections. 

3. LIFETIMES 

In this section we examine the probable coalescence 
times of relativistic double neutron stars, using a 
velocity-space contour map to illustrate dependence on 
the supernova kick velocity. As in the previous section, 
we avoid assumptions about specifics such as the distri- 
butions of stellar masses, kick velocities and orbital sep- 
arations, and present our results in units of the circular 
orbital speed and circular coalescence time (the merger 
time for an arbitrary system with zero eccentricity, de- 
noted as r c i rc ). These generic units have been chosen for 
their usefulness in comparing binaries in a scale-invariant 
manner. 

In producing the aforementioned contour map of co- 
alescence times, we first generated a family of arbi- 
trary double neutron star systems over all possible post- 
supernova relative orbital velocities. Computationally, 
this was achieved by uniformly populating the "velocity- 
space" with grid-points, where each point represents a 
relative orbital velocity. These velocities were then used 
to calculate the corresponding semi-major axes and ec- 
centricities, using Eq. and EJ Both of these factors 
will greatly influence the rapidity of each system's or- 
bital decay, with close, highly eccentric binaries being 
much more volatile than wider, more circular ones. To 
compute numerical estimates of just how rapidly each 
system will coalesce, each grid point's initial orbital pa- 
rameters were evolved forward in time by numerically 
integrating th e coupled differential equations derived by 
iPetersI l)1964f) for describing orbital decay by relativistic 
gravitational radiation. A Runge-Kutta 4th order algo- 
rithm was implemented with an adaptive step size, de- 
signed to integrate the equations forward to the precise 
time when the semi-major axis of the system drops to 
zero, thus determining the coalescence time. After it- 
erating through all grid-points, the 2D coalescence time 
data was used to produce the contour map shown in Fig. 

El 

The most significant feature of Fig. [21 is the large cen- 
tral region within which coalescence times rapidly drop 
from 0.1 T c i rc (in units of the coalescence time of a cir- 
cular system) to 10^ 5 r C j rc and even lower. To put this 
in perspective, for the orbital parameters of PSR J0737- 
3039A/B, 10~ 5 r C i rc corresponds to a coalescence time of 
less than 1,000 years. Fig. |3|thus demonstrates that a 
double neutron star system which receives a supernova 
kick of the order of its circular speed has a very high like- 
lihood of coalescing within an incredibly short timespan. 

The cumulative distribution of coalescence times is 
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Radial velocity component 

Fig. 3. — Contour map of coalescence times of binaries (in units 
of T c i rc ) as a function of orbital velocity immediately after the sec- 
ond supernova. Only prograde orbits are shown, but the diagram 
is symmetric about the x axis as for Fig 1. The time units are rela- 
tive to the coalescence time due to gravitational wave emission for 
a system with a circular orbit. If neutron stars receive a randomly 
oriented kick after the supernova explosion that produces them, 
it is possible to populate the entire phase space shown. Binaries 
with little angular momentum after the explosion are seen here to 
coalesce in remarkably short times relative to those in near-circular 
orbits. The consequence of this is that such systems are extremely 
difficult to observe compared to the systems of small eccentricity. 
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Fig. 4. — The fraction of systems that remain uncoalesced as 
a function of time, when considering a family of arbitrary binaries 
with all possible post-supernova orbital velocities. The solid curve 
here illustrates the case where the post-SN velocity is constrained 
to the pre-SN orbital plane, and the dashed line represents the 
general case. The unit used for time (T c i rc ) is the coalescence time 
of a system which remains in a circular orbit immediately after the 
supernova. For a sense of scale, consider that 1 r c ; rc is equivalent 
to 88 Myr when scaled to the parameters of PSR J0737-3039A/B. 



provided in Fig. 0] which shows the uncoalesced per- 
centage of the representative population as a function 
of time. This Figure dramatically illustrates the high 
probability of rapidly coalescing double neutron star sys- 
tems, with up to ~24% of all systems coalescing in less 
than 10~ 4 T C i rc (equivalent to 8,500 years when scaled 
to the parameters of PSR J0737-3039). After just 0.1 
Tcirc about half of all systems have coalesced. This is an 
intriguing result, implying that the coalescence rate of 
double neutron stars is underestimated. These rapidly 
coalescing systems would be practically undetectable by 
modern pulsar surveys, since they would not live long 
enough to give sufficient opportunity for detection. They 
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Fig. 5. — The evolution of observed eccentricities as a function 
of time, for the binary population as in Fig. [2] The units are again 
in terms of the time it would take a circular system to decay due to 
gravitational radiation (r c i rc ). After just 0.01 r c i rc the eccentricity 
distribution (dashed) has flattened significantly as all the short- 
period, highly eccentric systems have already coalesced. After one 
complete circular coalesence time, the population of surviving sys- 
tems (dotted) is much more uniform. 
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Fig. 6. — Orbital parameters of all currently known eccentric 
(e > 0.05) binary pulsar systems residing within the disk of our 
galaxy. We have excluded pulsars in globular clusters, where the 
formation histories are quite different. Relativistic systems (those 
which will coalesce within a Hubble time) are indicated by dot- 
ted circles, and non-coalescing systems are represented by stars. 
Note the comparatively low eccentricities of the relativistic bina- 
ries. P ulsar i nformation was obtained from the ATNF Pulsar Cat- 
alogue I Manchester ct al. 2005). 



would, however, be detectable by gravitational wave ob- 
servatories such as LIGO, which will come online in the 
near future. We might predict that at short orbital pe- 
riods there will be a deficit of highly eccentric systems 
because these are selected against due to their shorter 
lifetimes. We now investigate this in the following sec- 
tion. 

4. ECCENTRICITY DISTRIBUTION EVOLUTION 

Here we analyse the time evolution of the eccentricity 
distribution of a family of arbitrary double neutron star 
systems. The computation performed here is similar to 
that performed in the prior section - to begin, for every 
possible post-supernova relative orbital velocity (repre- 
sented computationally by uniformly spread grid-points 
in our "velocity-space" ) , we calculate an eccentricity and 
semi-major axis by applying Eq. This results in 

a standard representative family of double neutron star 
systems, whose orbital parameters we can then numeri- 
cally integrate forward in time us i ng th e coupled differen- 
tial equations derived bv iPetersI l)1964D . In this way, the 
eccentricity distribution of the whole family of possible 
double neutron star systems can be computed at birth, 
and at any time afterwards. Fig. [3] shows our result- 
ing computed eccentricity distributions at three stages 
in time. 

This Figure serves to illustrate the rapid coalescence 
of the vast majority of highly eccentric systems, which 
were produced in the greatest numbers when all possible 
bound orbital velocities were represented in our sample. 
This produces a rapidly evolving eccentricity distribu- 
tion which changes on very short timescales (0.01 r c ; rc ), 
and generally shifts towards lower eccentricities. While 
the initial post-supernova eccentricity distribution may 
be highly skewed, Fig. demonstrates that it doesn't 
retain its shape for very long, quickly flattening out to 
a more uniform distribution, from which the observed 
eccentricity distribution of double neutron star systems 



could be more plausibly derived. 

5. OBSERVED POPULATION 

Here we will examine the observed population of dou- 
ble neutron star binaries, particularly in terms of the 
eccentricity distribution. Fig. shows the semi-major 
axes and eccentricities of all currently known eccentric 
double neutron star systems residing within the disk of 
our galaxy. 

The majority of relativistic binaries tend to have low 
eccentricities, which seems fortuitous considering the 
tiny size of the low-eccentricity "bullseye" region in the 
velocity-space as shown in Fig. ^ To achieve low ec- 
centricities, these relativistic binary systems would have 
had to receive consistently precise supernova kicks. Any 
small deviation in kick direction or magnitude would re- 
sult in higher eccentricities. 

Our analysis of uniformly distributed pulsar kicks 
finds that surviving binaries should create mostly highly- 
eccentric systems. We contend that supernova kick veloc- 
ities can easily be large enough in magnitude to kick a bi- 
nary out of its small low-eccentricity "bullseye" and into 
the dominant, high-eccentricity region of velocity-space 
(see Fig. [y. These large supernova kicks can be inferred 
from the abundance of single pulsars roaming the galaxy 
at high velocities (mean of order 400-500 km s _1 ), most 
likely accelerated by supernova kicks. With regards to 
our simple kick model, it is almost certain that the real 
kick velocity distribution would not uniformly populate 
the available phase space after a supernova explosion, as 
the kicks in close binaries would have to range from small 
velocities to higher than 1000 km s _1 with equal prob- 
ability. Nevertheless, trying to guess the exact shape 
and magnitude of the pulsar kick velocity distribution 
leads to a plethora of models, and we resist this temp- 
tation. Our ignorance of the kick velocity distribution 
does not invalidate our results, unless pulsars can never 
receive a sufficiently large kick, of magnitude similar to 



() 



Chaurasia & Bailes 



that of the original pre-supernova orbital velocity. We 
are inclined to assume that such large kicks can exist, 
however, as observed pulsar velocities suggest kicks may 
be impa rted with magnitud es up to and exceeding 1000 
km s" 1 l|Cordes et al.lll993jh 

We propose that the observed population of relativistic 
double neutron stars exhibits low eccentricities because 
systems with higher eccentricities have coalesced before 
they could be detected. An average pulsar is detectable 
by modern pulsar surveys for approximately 10 Myr - 
however, we have shown that binary systems have an 
appreciable likelihood of being kicked into coalescence 
within 1 Myr or even less. Not only will these systems 
coalesce rapidly, they will also have orbits so tight that 
the effects of "Doppler smearing" are too strong for cur- 
rent search algorithms to detect them in some surveys. In 
summary, the low eccentricities of the relativistic double 
neutron star systems suggest that a strong selection ef- 
fect is in operation, and we claim that accelerated decay 
d ue to gravitatio n al rad iation is the cause. 

Faulkne r et alJ (J2005) suggest that there is a correla- 
tion reinforced by the discovery of PSR J1756-2251 be- 
tween orbital eccentricity and spin period of recycled pul- 
sars. Our results suggest that many highly-eccentric sys- 
tems resulting from close binaries are never seen by pul- 
sar surveys, and that long-period binaries favour high ec- 
centricities. Close binaries transfer matter when the evo- 
lutionary timescale of the donor is still long. We therefore 
expect that short-period progenitors will be spun up to 
shorter periods than the longer-orbital period systems. A 
natural consequence of this is that the short-period high 
e systems are selected against, that very few long-period 
low e systems are born, and therefore a correlation is 
observed. 

6. DISCUSSION: IMPLICATIONS FOR THE MERGER RATE 
AND "DOUBLE SUPERNOVAE" 

We have seen that the binary pulsar population is miss- 
ing highly-eccentric binary pulsars in relativistic orbits, 
which we attribute to their short observable lifetimes. 
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The most reliable estimates of the merger rate of neu- 
tron star binaries come from extrapolation of the ob- 
served population, which ignores these missing binaries 
that we expect to be formed in large numbers before 
dying young. Our calculations suggest that factors of 
several are required to correct the merger rate for these 
missing binaries. 

We have also shown that there is a non-negligible po- 
tential for coalescence times of the order of months if the 
pulsar kick is of the same order as the pre-supernova 
orbital velocity. This may mean that even relatively 
wide pre-supernova orbits can be driven to coalesce in 
remarkably short times by receiving an appropriately- 
oriented retrograde kick. This phenomenon could reveal 
itself in the form of a supernova followed closely by a sec- 
ond similarly-scaled catastrophic event (the coalescence) , 
where both would be observed at the same position a few 
months or years apart - what might be described as a 
"double supernova". How soon after the original burst 
of star formation in the Universe might we expect to see 
such an event? Population III stars - the ancient stars 
which were the first to form in the history of the Universe 
- may have often occurred in binaries where both stars 
underwent supernovae, subsequently taking part in an 
ultra-rapid coalescence which could be detected through 
gravitational radiation emissions (the only known way of 
observing these objects, otherwise hidden by the Epoch 
of Reionisation). This idea has been alluded to before 
l)Belczvnski et al.ll2004j) . however our findings emphasise 
the possibility that such a coalescence could occur very 
quickly after initial star formation, allowing future grav- 
itational wave detectors to directly study events that oc- 
curred within the first few Myr of star formation in the 
Universe. 
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